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The angular momentum projected Generator Coordinate Method using 
the quadrupole moment as collective coordinate and the Gogny force as the 
\ effective interaction is used to describe the properties of the ground state 

• and low-lying excited states of the neutron rich Hght nuclei 32, 34, 36 

' is found that the ground state of the nucleus ^'^Si is spherical. However, 

. this is not only due to the = 20 shell closure as the ground state of 

'^'^Si contains a significant amount of the intruder f-^ j2 neutron orbital. On 
the other hand, rather good agreement with experimental data for many 



o 



■ observables is obtained. 

7j ■ PACS numbers: 21.60.Jz, 21.60.-n, 21.10.Re, 21.10.Ky, 21.10.Dr, 27.30. +t 
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1. Introduction 



In nuclear physics the mean field approximation is always the first step 
to understand the properties of the ground and lowest-lying excited states. 
\ The mean field approximation provides the concept of magic numbers as 

well as the concept of spontaneous symmetry breaking. For nuclei with 
proton and/or neutron numbers close to the magic ones one expect symme- 
try conserving (i.e. non superconducting and spherical) ground states. On 
the other hand, for nuclei away from the magic configurations one expect 
strong symmetry breaking and the appearance of deformed ground states 
that generate bands (as the rotational bands). 

The experimental studies of light nuclei away from the stability line 
N = Z seem to imply that for those (usually neutron rich) nuclei some of the 
properties associated to magic numbers are not preserved. The most striking 
example is the experimental evidence towards the existence of quadrupole 
deformed ground states in the neutron-rich nuclei around the magic number 
N = 20. In addition, the extra binding energy coming from deformation 
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can help to extend thereby the neutron drip hne in this region far beyond 
what could be expected from spherical ground states. Among the variety 
of available experimental data, the most convincing evidence for a deformed 
ground state in the region around = 20 is found in the ^'^Mg nucleus 
where both the excitation energy of the lowest lying 2"*" state [1] and the 
B{E2, 0'^ — > 2+) transition probability [2] have been measured. Both quan- 
tities are fairly compatible with the expectations for a rotational state. At 
the mean field level, the ground state of ^^Mg is spherical. However, when 
the zero point rotational energy correction (ZPRE) is considered, the energy 
landscape as a function of the quadrupole moment changes dramatically and 
^"^Mg becomes deformed [3, 4, 5, 6, 7, 8, 9]. 

A more careful analysis of the energy landscape including the ZPRE cor- 
rection reveals that, in fact, there are two coexistent configurations (prolate 
and oblate) with comparable energy indicating thereby that configuration 
mixing of states with different quadrupole intrinsic deformation has to be 
considered. Therefore, an Angular Momentum Projected Generator Co- 
ordinate Method (AMPGCM) calculation with the quadrupole moment as 
collective coordinate is called for. We have applied this method in Ref. [10] 
to the study of the nuclei ^^~^'^Mg . The Gogny force [11] (with the DIS 
parameterization [12]) has been used in the calculations. We have obtained 
prolate ground states for •^'^^^'^Mg indicating that the = 20 shell closure 
is not preserved for the Mg isotopes. Moreover, a good agreement with 
the experimental data for the 2'^ excitation energies and B{E2) transition 
probabilities was obtained. The purpose of this paper is to extend those 
calculations to the study of the nuclei ^^"•^^S'i which are obtained by adding 
two protons to the ^^'^"^Mg nuclei. These nuclei have recently being studied 
experimentally [13] and their ground states seem to be spherical and there- 
fore preserve the N = 20 shell closure. Therefore, our intention is to check 
whether our method and interaction is not only able to predict a deformed 
ground state in ^^Mg (Z=12) but also to predict an spherical one in ^Si 



The angular momentum projected Generator Coordinate Method (AMP- 
GCM) with the mass quadrupole moment as generating coordinate is used 
as the theoretical framework. As we restrict ourselves to axially symmetric 
configurations, we use the following ansatz for the K = wave functions of 
the system 



(Z=14). 



2. Theoretical framework 




(1) 
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In this expression \^{q2o)) is the set of axially symmetric (i.e. = 0) 
Hartree-Fock-BogoHubov (HFB) wave functions generated by constraining 
the mass quadrupole moment to the desired values ^20 = {'■p{'l2o)\ — 
l/2(.x^ + W{(l2o))- The quasiparticle operators associated to the intrinsic 
wave functions \^{q2o)) have been expanded in a Harmonic Oscillator (HO) 
basis containing 10 major shells and with equal oscillator lengths to make 
the basis closed under rotations [14]. The operator 

^00 = ^^^^ / dnd'Me-^^^^e-'P^ye-^^^' (2) 

is the usual angular momentum projector with the K = Q restriction [15] 
and fi{q2o) are the "collective wave functions" solution of the Hill- Wheeler 
(HW) equation 

j dq2on\q20,q2o)fai'l2o) =-^^/ dq20-^\q2Q,q2o)f!j{<l2o)- (3) 

In the equation above we have introduced the projected norm 

M\q20. qho) = {^{Q2o)\ Pio Hq2o)) (4) 
and the projected hamiltonian kernel 

n'{q20,q2o) = {v{q2o)\HP^o Hqho)) ■ (5) 

As the generating states Pjo \ ^{q2o)) are not orthogonal, the "collective am- 
plitudes" fl{q2o) cannot be easily interpreted. Instead, one usually introduce 
[16] the "collective" amplitudes 

which are orthonormal / dq2o{gi)*{q2o)g^'iq2o) = 5i,i'^<j,(j' and therefore 
their module squared has the meaning of a probability. 

The B{E2) transition probabilities are computed using the AMPGCM 
wave functions as 

B{E2,I, ^ If) = j dq^odqLfaj*{qL){IflL II Q2 II hq^o)faM\o) \ 

(7) 

As we are using the full configuration space no effective charges are needed. 
Further details on the computational procedure can be found in Ref [10]. 
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3. Discussion of the results 



In Figure 1 we show the HFB potential energy surfaces (PES) for the 
three nuclei considered (dashed line) along with the projected energies = 
^^(92, 92) /AA^ (92, 92) for I = 0, 2, 4, 6 and 8 (full lines). The HFB PES show 
a very pronounce spherical minimum in ^"^Si, a shallow spherical minimum 
in ^^Si and, finally, a very swallow oblate minimum in These results 

are rather different from the ones obtained for the corresponding Mg iso- 
topes with the same neutron number [10]. The reason for such differences 
is that the Si isotopes have the proton shell completely occupied and 
therefore, the up-slopping K = 5/2 level of the proton orbital prevents 
the appearance of prolate deformation. At the mean field level we can say 
that the Si isotopes around = 20 retain the closed shell properties associ- 
ated to the = 20 shell closure. However, as it has already been recognized, 
the effect of the restoration of the rotational symmetry in the HFB PES can 
be very important and it can substantially modify the conclusions extracted 
from the mean field results. The angular momentum projected (AMP) PES 
of Figure 1 show in the three nuclei considered two minima, one prolate and 
the other oblate, for all the angular momentum considered. However, the 
prolate and oblate minima lie rather close in energy in almost all the cases 
and therefore, we expect strong mixing of the two configurations. Taking 
into account this fact and also that the amount of mixing is not only deter- 
mined by the AMP PES but also by the "collective inertia" the amount of 
mixing can only be disentangled when a configuration mixing calculation is 
carried out. 




Fig. 1. The HFB potential energy surface (PES) for the three nuclei considered 
(dashed line) together with the projected energies for / = 0, . . . , 8 (full lines). 
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We have also carried out configuration mixing (in the framework of the 
GCM) calculations taking into account angular momentum projected wave 
functions and the quadrupole constrained HFB wave functions as generating 
intrinsic states. The solution of the HW equation provides us with energies 
E"^ and collective wave functions /J(g2o)- As mentioned in section 2 the 
amplitudes are not well suited for a physical interpretation and there- 
fore one considers the collective amplitudes gi{q2o) instead. In terms of 
these amplitudes one can defined an "averaged" intrinsic quadrupole mo- 



ment (gf2o)i = Idq2oq20 giiq 



220 J 



for each of the AMPGCM solutions. 
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Fig. 2. The HFB potential energy surface (PES) for the three nuclei considered 
(dashed line) together with the projected energies for I = 0,2,4,6 and 8 (full 
lines). 

In Figure 2 we show, along with the angular momentum projected PES 
(plotted again to guide the eye), the energies E^. solution of the HW equation 
for (7 = 1 and 2 and angular momenta / = 0,2,4,6, and 8. Each energy 
has been placed at a value of the quadrupole moment corresponding to the 
"averaged" intrinsic quadrupole moment {q2o)i- From figure 2 we see that in 
^^Si the two 0"*" states shown are spherical, the 2f and 4f are oblate whereas 
the 2 J, 4^, 6t and 8^ are prolate deformed states. The 6^ state is almost 
degenerate in energy with the nearly spherical Qf state and this could explain 
the anomalous quadrupole moment of the former state. A nearly identical 
pattern is also seeing in the nucleus ''''^Si except for the position of the 0^ 
state that lies below the 2'^ . In the nucleus ^^Si the ground state is spherical 
and a prolate deformed 0^ state is obtained. The 2+ and 4"*" states are all of 
them almost spherical whereas the 6'^ and 8]*" are prolate and the 6^ and 83" 
states are oblate. We observe that in the three nuclei considered the ground 
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state is spherical due to the configuration mixing between the prolate and 
oblate minima. The rest of the spectra shows rather unclear patterns that 
can only be elucidated by looking at the "collective amplitudes" g^{q2o)- 
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Fig. 3. The collective amplitudes 15^(920)^ (thick lines) for tr = 1 (full) and 2 
(dashed) and spin values of J = 0?i, . . . , 8?i for the nuclei ^'^Si, ^^Si and ^^Si. The 
projected energy curve for each spin is also plotted (thin Hue). The y-axis scales are 
in energy units and span an energy interval of 13 MeV. The collective amplitudes 
have also been plotted against the energy scale after a proper scaling and shifting, 
that is, the quantity E^^ + 15 x 15^(920)^ is the one actually plotted. 
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In figure 3 the collective wave functions squared gi{q2o) for the two 
lowest solutions a = 1 and 2 obtained in the AMP-GCM calculations are 
depicted. We also show in each panel the projected energy for the corre- 
sponding spin. We observe that the 0^ ground state wave functions of the 
three nuclei contain significant admixtures of both the prolate and oblate 
configurations. Looking into the single particle diagram of ^'^Si (not shown 
here) it is observed that the neutron /7/2 shell, that is empty at q2o = 
and 5.78 MeV higher than the crosses the Fermi surface at the defor- 
mations q2o = —0.56 in the oblate side (the K = 7/2 level) and q2o = 0.56 
in the prolate one (the K = 1/2 level). Therefore, the intrinsic states with 
q2o > 0.56 and q2o < —0.56 have (at least) two particles in the intruder fj/2 
orbital. As the ground state collective amplitude in ^"^Si (and in the other 
two nuclei) is not negligible in the intervals q2o > 0.56 and q2o < —0.56 
we conclude that the ground state of the three nuclei contain significant 
amounts of the intruder fi/2 neutron configuration in spite of being spheri- 
cal. To make the argument more quantitative we have computed the quan- 

r 2 

tity 1 — /_o^5 dq2o gi^^{q2o) that gives us an idea of the percentage of the 
"collective" amplitude outside the ^20 interval between -0.5 and 0.5 6 (i.e. 
the interval where the intrinsic states contain the intruder configura- 
tion). This quantity is roughly 0.35 for o" = 1 in the three nuclei considered 
and it corresponds roughly to 0.7 particles in the neutron 77/2 sub-shell. 
Therefore, we conclude that the ground state of ^Si does not correspond 
to a pure configuration with A'' = 20 acting as a closed shell. Similar con- 
clusions have been obtained in shell model calculations [17] and also in the 
Shell Model Monte Carlo calculations of [18]. 

The collective amplitude of the 02^ excited state of ^Si is rather unusual, 
being concentrated around ^20 = 0. This unusual behavior probably is the 
responsible for the low excitation energy of this state (see below). In the 
^^Si and ^Si nuclei and 1 = 2 and 4 the collective amplitudes are well 
located inside the corresponding prolate and oblate minima. For I = 6 we 
have accidentally coexisting prolate and oblate minima and for I = 8 the 
lowest minima becomes the prolate one. In the nucleus ^^Si and 1 = 2 and 
4 we have shape coexistence between the prolate and oblate minima and 
therefore the collective amplitudes are equally distributed among them; the 
net result being that these states are, on the average, spherical. For 1 = 6 
and 8 the prolate minimum becomes dominant and the collective amplitudes 
are localized around the minima. 

Now, we want to compare our results with the available experimental 
data. The first piece of information we can compare to is the two neutron 
and two proton separation energies. The two neutron separation energies 
{S(2n)) for ^"^Si and ^^Si are 11.83 MeV and 9.6 MeV respectively. These 
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numbers have to be compared with the experimental data [19]: 12.018 MeV 
and 8.587 MeV, respectively. We observe a reasonable good agreement with 
experiment for both nuclei. Concerning the two proton separation energies 
{S{2p)) we have used the previous results [10] for the magnesium isotopes 
to obtain 28.70 MeV, 32.72 MeV and 36.26 MeV for '•^'^-'■^^ Si, respectively. 
The corresponding experimental values [19] are 29.777 MeV, 33.737 MeV 
and 35.430 MeV. The agreement between theory and experiment is again 
rather reasonable. 





Energies (MeV) 


Exp. 


S(£;2)e^/m'* 


Exp. 




a 


b 


c 


a 


d 


e 


f 


d 


''''Si 


2.34 


3.60 


5.95 


1.941 


82.18 


0.01 


28.00 


113 ± 33 


'■'''Si 


2.92 


2.52 


5.62 


3.327 


108.99 


50.41 


76.92 


85±33 


^""Si 


2.63 


3.24 


4.32 


1.399 


211.55 


65.44 


48.81 


193 ± 59 



Table 1. Calculated and experimental results for excitation energies (in MeV) and 
B(S2,0+ — > 2+^) transition probabilities (in e^/m^). The columns marked a, b 
and c correspond to 0^ — 2^,0^ — 0^ and 0^ — 2^ respectively. The ones marked 
d,e and f correspond to Of — 2f , O2 — '^i and 0^ — 2 J respectively. 

In table 1 the energy splittings between different states and the B{E2) 
transition probabilities among some of them are compared with the available 
experimental data. Concerning the B{E2,{)'f 2J^) transition probabilities 
we find a quite good agreement with the experiment. In fact, our calculation 
reproduces the increase of the B{E2) values in going from ^'^Si to ^^Si. The 
other transition probabilities are much smaller than the 0]*" — > 2^ ones ex- 
cept in ^Si. On the other hand, the theoretical values for the 2^ excitation 
energies agree reasonably well with the experimental data for ^'Si and ^'^Si 
but shows a discrepancy of 1.2 MeV in the case of '^^Si. This discrepancy can 
be attributed to many sources, like not taking into account triaxial intrinsic 
configurations or not dealing with pairing correlations in a "beyond mean 
field" framework. Another source for the discrepancy could be related to 
the fact that ours is a calculation of the Projection After Variation (PAV) 
type instead of the more complete projection before variation (PBV). Usu- 
ally, the PBV method yields to rotational bands with moments of inertia 
larger than the PAV ones [20, 21]. In Ref [10] we estimated the effect of 
considering PBV in calculations similar to the ones performed here but for 
some deformed magnesium isotopes. The conclusion was that the moment 
of inertia gets enhanced by a factor 1.4 and therefore the excitation energies 
have to be quenched by a factor 0.7. Unfortunately, the method used in [10] 
for the mentioned estimation can not be applied here as we are dealing with 
spherical ground states as well as with nearly spherical excited states. 
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On the other hand, shell model calculations [13, 17, 22] predict a second 
O2" excited state in ^Si with an excitation energy (2 MeV in [13], 3 MeV 
in [17] and 2.6 MeV in [22]) which is below the excitation energy of the 2^ 
state. In the experimental work of Ref [13] such excited state was not ob- 
served but it was argued that with the given experimental setup it was quite 
difficult to detect it. However, in the experiment reported in [22] the men- 
tioned excited 0"*" state was observed at a tentative excitation energy of 2.1 
MeV. In our calculations we also observe a second 0^ state in ^'^Si (spherical 
in character) which lies at an excitation energy of 2.52 MeV. Our predicted 
B{E2) transition probability 2^ is 50.41e^/m^ which is a factor four 

smaller than the one obtained in the shell model calculations of [13, 17] and 
a factor of six smaller than the predicted value of [22]. In [22] an exper- 
imental indirect estimate of the B{E2) transition probability of 444(210) 
e^fm^ is given. At this point it has to be recalled that the B{E2) tran- 
sition probabilities are more sensitive to the collective wave functions than 
the excitation energies. Therefore, a relatively small change in the collective 
wave functions can alter drastically the transition probabilities. In our case 
the collective amplitude of the 0^ state (see Fig. 3) is concentrated around 
the spherical intrinsic configuration (result which is no consistent with the 
shell model expectations of [17, 22]). This spherical intrinsic configuration 
has zero pairing correlations in our HFB calculations and therefore it can be 
expected that the inclusion of pairing correlations in a "beyond mean field" 
framework can modify the properties of the quadrupole dynamics around 
this point. To elucidate this and other deficiencies a simultaneous projec- 
tion on angular momentum and particle number (like to one performed in 
[8] with the Skyrme interaction) is called for. 



4. Conclusions 

In conclusion, we have performed angular momentum projected Gener- 
ator Coordinate Method calculations with the Gogny interaction DIS and 
the mass quadrupole moment as generating coordinate in order to describe 
quadrupole collectivity in the nuclei ^'^Si, ^Si and ^^Si. We obtain a spher- 
ical ground state in ^^Si (in opposition to the nucleus ^"^Mg) but it contains 
significant admixtures of the intruder configuration with two neutrons in 
the /■j/2 orbital. We also find in ^Si that the first excited state is a 0+ 
with an excitation energy of 2.5 MeV. In all the nuclei considered the excita- 
tion energy of the lowest 2~^ state is rather well reproduced. Moreover, the 
B{E2) transition probability from the ground state to the 2+ state is also 
well reproduced in the three cases. Other quantities like two neutron and 
two proton separation energies also compare well with experimental data. 
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